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We have examined modifications to the nuclear weak charge due to small differences 
between the spatial distributions of neutrons and protons in the Cs nucleus. We derive 
0^ ' approximate formulae to estimate the value and uncertainty of this modification based 

only on nuclear rms neutron and proton radii. Present uncertainties in neutron distribu- 
tions in Cs are difficult to quantify, but we conclude that they should not be neglected 
when using atomic parity nonconservation experiments as a means to test the Standard 
^ ' Model. 
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' Recent measurements of transition polarizabilities |^ , coupled with previous measurements of parity 
, nonconservation (PNC) in atomic cesium have significantly reduced uncertainties associated with 
the extraction of Q^,, the radiatively corrected weak charge of the Cs nucleus. The latest result 
QJ^P* = — 72.06(28)cxpt (34)atomic theory IS in mild disagreement, at the 2.5cr level, with the Standard 
' Model prediction of Q^'^"'^- — — 73.20(13)thoory |§1 The experimental number requires input from atomic 
\Q . theory calculations which include effects of normalization of the relevant axial electron transition 
matrix element in the vicinity of the nucleus. The finite nuclear size is incorporated by including pN{r), 
the spatial nuclear distribution, in the matrix elements. One possible contribution to Qu, which has been 
left out of the quoted numbers is the modification of the extracted weak charge due to the difference 
. between neutron and proton spatial distributions in this nucleus with relatively large neutron excess. 

The effect of the neutron distribution differing from the proton distribution in a nucleus has been 
explicitly considered in the atomic theory calculations, [Q and was dismissed because the estimated 
size was extremely small compared to existing uncertainties at the time. Other authors |^,^ have also 
^ , derived and discussed this contribution further. In the case of Cs, all authors agree the effect is quite 
^ ' small. However, with the significant reduction in errors in recent atomic PNC measurements, the effect 
should no longer be neglected. As we argue below, the additional uncertainties in extracting from the 
data arising from neutron-proton distribution differences are slightly below the uncertainties arising from 
, atomic theory calculations or current experimental error bars, but are comparable to Standard Model 
radiative correction uncertainties. 

In this note, we attempt to quantify the additive contribution and uncertainties to the nuclear weak 
charge, Quj, arising from the relatively poorly known spatial distribution of neutrons in the nucleus, 
Pn{f)- We briefly summarize some relevant nuclear structure issues, both theoretical and experimental. 
We also briefly discuss methods that could improve this knowledge. We present results of our numerical 
calculations of Qw arising from various p„ distributions, and present approximate methods which show 
what effect differing nuclear structure model predictions would have on precision Standard Model tests. 

At tree level in the Standard Model, the nuclear weak charge is Q^'^°'^ = (1 — 4sin^ Ow)Z — N, with 
N and Z the neutron and proton number, and sin^ 6w the weak mixing angle. Standard Model radiative 
corrections modify this formula slightly. ||] The effect of flnite nuclear extent is to modify N and Z to 
QnN and qpZ respectively, [|j where 

qn{p) = J f{r)p„(p){r)d^r. (1) 



> 



O 

OS 



I 



*e-mail: Steven.Pollock@Colorado.edu 
t e-mail : welliver @ lucky. Colorado . edu 



1 



Here /(r) is a folding function determined from the radial dependence of the electron axial transition 
matrix element inside the nucleus, and the neutron (proton) spatial distribution Pn(p) is normalized to 
unity. It is common to characterize the neutron distribution by its rms value, i?„, since it can easily 
be shown that the weak charge is most sensitive to this moment. To the extent that p„ and pp are the 
same, the overall nuclear size effect can be completely factored out. This has explicitly been done in the 
experimental extraction of Qw . The slight difference between g„ and qp has the effect of modifying the 
effective weak charge: 

Q^ = QSt.Mod^^g„-p^ (2) 

where Q 

AQl-P = Nil - q,Jqp). (3) 

A naive calculation, helpful for quick estimates of the effect of different possible neutron distributions 
on AQJ^~P, can be made by assuming a uniform nuclear charge distribution (zero-temperature Fermi gas), 
and then parameterizing the neutron distribution solely by its value of i?„. In this approximation, one 
solves the Dirac equation for the electron axial matrix elements, /(r), near the origin by expanding in 
powers of a (the fine structure constant). Finally, we can assume Rn ~ Rp, characterizing the difference 
by a single small parameter, (R^/Rp) = 1 + e. In this case, we find |^ 

qp^ l-(Za)2(.26), (4) 
g„«l-(Za)2(.26 + .221e), (5) 
AQl-P « N{Za)^.221e)/qp. (6) 

Eq. ^ shows the rough dependence of the correction to the weak charge on the difference between 
neutron and proton distributions, characterized by e. Results of this naive calculation are shown as the 
solid line in Fig. |l|. The slope of the line demonstrates the sensitivity of the uncertainty in AQ^~p to the 
uncertainty in rms neutron radius. The range in e of ±0.1 corresponds to a SRn/Rn of about 5%, which 
we argue below might be a reasonable estimate of the uncertainty in neutron rms radius. We do not 
need to rely on these approximations; we have solved the Dirac Equation numerically for s- and p- state 
electron wave functions given the experimental charge distribution of Cs, evaluated /(r) numerically, 
and thus calculated g„, qp, and AQ^~p given various model predictions for the neutron distribution. 
The diamonds in Fig. |^ correspond to the full calculation assuming neutron distributions with the same 
shape as the proton distribution, scaled to give the particular values of e. The above approximations 
prove to be accurate, although the resulting uncertainty in AQ^~p is marginally underestimated by only 
including the uncertainty in The additional effects of neutron distribution shape variations could 
slightly increase the uncertainty in the nuclear contribution to the weak charge, as demonstrated by the 
error bars on the diamonds in Fig. |l|. These error bars arise by assuming a 2 parameter Fermi fit for the 
neutron distribution, Pn{r) — 1/(1 + e['^''~^-'/^l), and allowing the "skin thickness" parameter z„ to vary by 
±0.1, keeping e fixed. Such a variation is comparable to the difference z„ — Zp in various nuclear models 
[Q. (A more detailed analysis of the effect of neutron shape on will be presented elsewhere 0.) 
From Fig 1, it is clear that the uncertainty in the radius i?„ dominates the uncertainty in AQ""^. 

The effect of AQ^~p was understood and estimated in the atomic structure calculations Q by first 
assuming p„(r) = Pp{r) (so qn = qp and AQ^~p = 0) and then recalculating with a theoretical param- 
eterization of the neutron density, jsj The resulting AQ^~p sa 0.06 was extremely small, amounting to 
about 0.08% of the total weak charge, and was thereafter ignored. This neutron density was obtained by 
scaling a variational extended spherical Thomas-Fermi calculation using an effective parameterization of 



'^There will be additional small multiplicative corrections to AQ'^~'' arising from Standard Model radiative 
corrections, as well as additive corrections arising from e.g. internal structure of the nucleon, but these can be 
safely neglected since AQ^~^ is itself so small. 
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the nuclear Lagrangian ( "Skyrme SkM*" ) , which happened to yield a nuclear neutron rms radius which 
differed by only 0.9% from the proton rms radius. Eq. |^ confirms the size of this shift, given only the rms 
neutron and proton radii. However, the assumed Pnir) distribution may not be an accurate representa- 
tion of the correct neutron distribution. There exists both theoretical and experimental evidence that 
Rn might differ from Rp by significantly more than 0.9%, and thus from Eq. |, AQ^,"^ 

may be similarly 

underestimated. 

In a more recent theoretical analysis, Chen and Vogel 0| considered two more sophisticated nuclear 
structure models. Both models involved a Skyrme parameterized nuclear Lagrangian, |^ computed in 
the spherical Hartree-Fock (HF) approximation. Such models are quite successful in predicting a wide 
variety of nuclear observables, including charge distributions, binding energies, bulk properties, etc. 
These two models (SkM* and Skill) yielded Rn/Rp values of 1.022 and 1.016 respectively. Using the 
average of these values in Eq. ^ we obtain AQJ^'P = +.11, double the estimate of Ref. [Q. Using spatial 
distributions of neutron and proton densities from even more recent nuclear structure models , we have 
calculated the nuclear correction directly, rather than using the approximation of Eq. ^. Using spherical 
Skyrme SLy4 distributions, we find AQJ^^^ — +.14. Similarly, using (spherical) Gogny distributions, 
including blocking, we find AQl-P = +.11 (Eq. || for these two cases predicts +.15 and +.12 respectively). 
Relativistic potentials ||l0[| typically generate significantly larger neutron radii (see discussion below), and 
thus would predict larger AQJJ^^p, possibly by a factor of 2 or more, based on calculations in nearby nuclei, 
but no ^^'^Cs distributions for such models have been published to date. Note that if i?„ > Rp, then 
AQ^~P is positive. The central value of the most recent experiment 0] gives Q^, — —72.06, compared to 
QSt.Mod _ —73.20, so this nuclear correction is of the right sign to partially explain the small discrepancy. 
However, if one wanted to attribute the difference entirely to nuclear physics effects, one would require 
Rn = (1.18 ± .Q7)Rp (adding all atomic experimental and theoretical, and Standard Model theoretical 
errors in quadrature), which is significantly out of the range of any theoretical or experimental nuclear 
structure predictions. 

The fundamental question regarding nuclear structure remains — what uncertainty should be associ- 
ated with AQ^^P? Chen and Vogel argued that a reasonable uncertainty in their calculated neutron 
radius might be SR^^ w ±lfm^. According to Eq. ||, this corresponds to an uncertainty dAQ^~P = ±0.13. 
The estimate in ref. Q for the theoretical uncertainty in Q^, for a single isotope was slightly larger, 0.25% 
of Qu,, i.e. ±0.18. This is still quite small compared to the current atomic structure uncertainty (±0.34 
in Qw)i but is as large as the uncertainty in Q^, arising from uncertainties in Standard Model radiative 
corrections (see Fig. ^ . All of the models we have considered predict the charge radius in Cs within about 
1%, but the parameter fits used to determine the Skyrme potentials are based in part on observables, 
including charge radii, in nearby semi-magic even-even nuclei. There remain various possible sources of 
concern that a value of SR^ ±lfm^ may still be an underestimate. For example, ^'^^Cs is a deformed, 
odd-Z nucleus. Most nuclear structure calculations for large nuclei assume spherical symmetry with at 
least partially closed nuclear subshells. Pairing and blocking effects make calculations with odd N or Z 
less reliable, as evidenced by the failure of most Skyrme HF calculations to reproduce experimentally 
observed "even-odd" staggering of charge radius along isotope chains. |l^ In reference Q pairing effects 
were included, but deformation was included only in a semi-phenomenological manner. 

There exist other classes of nuclear structure models which give quite different predictions for neu- 
tron properties, for example, relativistic Hartree models based on a modified Walecka-model nuclear 
Lagrangian. |l0| These models have seen significant improvements in recent years, and may now be 
viewed as competitive with more established Skyrme models in terms of their predictive power over a 
wide variety of observables throughout the periodic table. In a recent paper comparing models, |pl| 
R^/Rp for ^^^Ba (the nearest even-even semi-magic nucleus above Cs) ranged from 1.03 in a Skyrme 
model to more than 1.08 in the relativistic models. The difference in predicted R^ between these two 
models alone exceeds 1 fm^. For ^^^Xe, R'^/Rp values vary from 1.04 to 1.09, with the predicted i?^ 
differing by well over 1 fm^. In another recent paper comparing models, jlTt the predictions for i?„ in 
^^''Sn (with a value of N/Z similar to "'^^'^Cs) varied by more than 2 fm^ between extreme models, a spread 
of over 8%. Again, these calculations are primarily for even-even nuclei; relativistic models have not yet 
been used to calculate self-consistently in the neighboring unpaired (odd Z) cases. This only adds to the 
uncertainty in the prediction of a model spread for the case of Cs. Based on these spreads, it appears 
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that current nuclear theory yields an uncertainty of at least 4 or 5% in i?„. 

The uncertainties in neutron distributions discussed so far arise from disagreements between model 
predictions. It is important to note that the neutron rms radius has never been directly measured in 
any isotope of Cs. Indeed, it is extremely difficult to measure i?„ in any nucleus - the most accurate 
measurements of charge radii come from electromagnetic interactions, which are dominated by the proton 
distribution. Elastic magnetic scattering is affected mostly by unpaired (valence) nucleons, which does 
not allow for a detailed or accurate measure of the bulk rms neutron radius. Data from strong interaction 
probes measure the "matter radius" , but are somewhat more sensitive to surface effects, and suffer from 
some poorly controlled systematic theoretical uncertainties arising from the models required in analyzing 
strong interaction observables. For example, there exist data from polarized proton elastic scattering on 
heavy nuclei, p^ ] The data are statistically of high quality, and are frequently viewed as an accurate 
experimental measure of i?„ in several heavy nuclei, including Sn and Pb. However, the systematic 
uncertainties in extracting i?„ , including choice of optical model and spurious variations in the result as a 
function of experimental beam energy, easily approach 5% or more. Other data, including pion or alpha 
scattering, suffer from similar uncertainties. The experimentally extracted average value from polarized 
proton scattering and pionic atoms for i?„ in ^ospj-, jjffgj. j^y around 3%. 

Even if strong interaction measurements can be argued to provide an accurate measure of the neutron 
rms radius, the weak interaction is sensitive to the spatial distribution of weak charge, which can not be 
exactly identified with neutrons or protons, but also includes effects of other nuclear degrees of freedom 
including e.g. meson exchange, and is more sensitive to non-surface density variations. A parity violating 
electron scattering experiment p^-[l6|| could directly measure the weak charge distribution, precisely 
what is needed for the interpretation of atomic PNC as a standard model test, and would be of clear 
value. Even if measured on another nucleus, the additional constraint on nuclear models should increase 
confidence in the predicted neutron distribution in Cs. As can be seen from Eq. ^, high precision atomic 
PNC measurements on significantly higher Z nuclei are more sensitive to the neutron distribution than 
in the case of Cs. Thus, a measurement of atomic PNC on extremely heavy nuclei might also be used as 
a measure of the neutron distribution, which in turn could be used to constrain the isovector parameters 
in the nuclear models, and thus increase the reliability of the predictions for Cs. 

To summarize, AQJJ,"^ is the deviation between the experimentally extracted weak charge and Stan- 
dard Model predictions due solely to differences in neutron and proton weak charge spatial distributions. 
The predicted value is small, typically of order 0.1, but with an uncertainty larger than the value it- 
self, arising mostly from uncertainties in i?„. This should be compared to the nominal value of the weak 
charge, Q^'^°'^ = —73.20(13). The effect of uncertain nuclear structure is thus comparable to the present 
uncertainties involved in the Standard Model prediction, and for i?„ > Rp is of the right sign to par- 
tially explain the experimental discrepancy in Cs. For these reasons, it should be included in any future 
atomic PNC tests of the Standard Model. With any significant further reduction in the uncertainties in 
atomic theory calculations, this nuclear contribution may eventually limit the level at which Standard 
Model tests can be performed with atomic PNC on Cs. To reliably reduce this uncertainty would require 
additional direct experimental input on neutron distributions, most likely from parity violating electron 
scattering at low momentum transfer, e.g. at a facility such as Jefferson Lab. 
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FIG. 1. Correction to the weak charge of Cs due to differences between neutron and proton spatial distributions 
as a function of e = R^/Rp — 1. The sohd curve is the prediction of Eq. ^ The diamonds are generated from a full 
numerical calculation with a realistic charge distribution, scaling pp to get p„. The error bars on the diamonds 
arise by varying the surface thickness of p„, keeping e fixed. The circle uses Eq. ^ assuming the neutron radius 
and uncertainty of ref. [Q| . The vertical error bar to the side of the plot shows just the uncertainty in the Standard 
Model prediction Q of Qw- 
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